Frequency dependence of induced spin polarization and spin current in quantum wells 
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Dynamic response of two-dimensional electron systems with spin-orbit interaction is studied theo- 
retically, on the basis of quantum kinetic equation taking into account elastic scattering of electrons. 
The spin polarization and spin current induced by the applied electric field are calculated for the 
whole class of electron systems described by p-linear spin-orbit Hamiltonians. The absence of 
non-equilibrium intrinsic static spin currents is confirmed for these systems with arbitrary (non- 
parabolic) electron energy spectrum. Relations between the spin polarization, spin current, and 
electric current are established. The general results are applied to the quantum wells grown in 
[001] and [110] crystallographic directions, with both Rashba and Dresselhaus types of spin-orbit 
coupling. It is shown that the existence of the fixed (momentum-independent) precession axes in 
[001]-grown wells with equal Rashba and Dresselhaus spin velocities or in symmetric [110]-grown 
wells leads to vanishing spin polarizability at arbitrary frequency uj of the applied electric field. 
This property is explained by the absence of Dyakonov-Perel-Kachorovskii spin relaxation for the 
spins polarized along these precession axes. As a result, a considerable frequency dispersion of spin 
polarization at very low u) in the vicinity of the fixed precession axes is predicted. Possible effects 
of extrinsic spin-orbit coupling on the obtained results are discussed. 

PACS numbers: 73.63.Hs, 72.25.Pn, 73.50.Mx 



I. INTRODUCTION 



The presence of spin-orbit interaction in solids provides 
a natural way of manipulating spin states of electrons by 
purely electrical means, without application of a mag- 
netic field. This property is a subject of interest for the 
novel and rapidly developing field of spintronics.^ Some 
important manifestations of spin manipulation are the 
generation of non-equilibrium spin polarization of elec- 
trons and excitation of spin currents by a driving elec- 
tric field which also leads to the usual charge current. 
Though the problem of field-induced spin polarization 
is an old one,^''^''*'^ it has been recently set at the fo- 
cus of attention. The main reason for this is the ap- 
pearance of experimental works^''''^'®'^*''^^ which demon- 
strate the spin polarization generated under the current 
flow both in bulk and two-dimensional (2D) semiconduc- 
tor layers with spin-orbit interaction. Another reason is 
the rapidly growing interest to the related phenomenon, 
the intrinsic spin-Hall effect, when the electric current 
in the presence of spin-orbit coupling leads to a non- 
equilibrium spin current in perpendicular direction (note 
that weak spin currents can exist even in equilibrium^^). 
After the theoretical proposal^^ of the intrinsic spin-Hall 
effect based on the Rashba model representing spin-orbit 
interaction for 2D electrons in quantum wells, it has 
been realized^^'i^-i'^'i^'i^'^O'^i that this effect is absent 
in the infinite 2D system in the static (zero-frequency) 
limit and exists at non-zero frequency or in finite-size 
samples. The absence of the static intrinsic spin-Hall ef- 
fect is not a general property, since it is related to the 
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specific form of the spin-orbit Hamiltonian. It has been 
shown20'22,23,24,25,26 ^-^g^^ g^j^y^ spiu-HaU effect ex- 
ists for more complicated models involving higher-order 
(cubic) terms in the momentum dependence of the spin- 
orbit Hamiltonian. The experimental observation of the 
spin-Hall effect for the 2D hole system^^ described by the 
spin-orbit Hamiltonian of this kind has confirmed this 
conclusion. The edge spin accumulation due to the spin- 
Hall effect has been also observed®'^^'^'* for 2D electron 
systems. Recently, it has been suggested^^''^" that devi- 
ation of electron band dispersion from the parabolic one 
can lead to non-zero intrinsic static spin-Hall effect. This 
point, however, remains controversial, since the analysis 
given in Ref. 17 predicts zero static spin-Hall effect even 
in this case. The calculations presented in this paper also 
show the absence of static spin currents for non-parabolic 
electron band dispersion. 

In contrast to the static regime, the frequency- 
dependent induced spin polarization and spin current 
have not been extensively studied. Theoretical calcula- 
tions of the frequency-dependent spin-Hall current based 
on the Rashba model have been done in Refs. 16, 18 
and 31 by using the methods of non-equilibrium Green's 
functions, Kubo-Greenwood linear response theory, and 
quantum kinetic equation, respectively, with the same 
result. The authors of Ref. 31 also calculated the fre- 
quency dependence of the induced spin polarization. A 
more complicated case of frequency-dependent response, 
when both Rashba and Dresselhaus (linear in momen- 
tum) terms are included in the spin-orbit Hamiltonian, 
has been studied in Refs. 32 and 33 in the dynamical 
(collisionless) regime. A comparative numerical study of 
the frequency dependence of the spin-Hall effect has been 
done in Ref. 23 for linear, cubic, and modified Rashba 
models. The resonances in frequency dependence of spin- 
Hall conductivity in magnetic field have been described in 
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Ref. 34 for the Rashba model. Nevertheless, a systematic 
investigation of the frequency-dependent problem of the 
induced spin polarization and spin current is still missing. 

A step towards systematic description of the frequency- 
dependent spin response is undertaken in this paper by 
considering the important class of spin-orbit Hamiltoni- 
ans: 



hp = hilp 



(1) 



where & is the vector of Pauli matrices, p = {px,Py) is 
the 2D momentum of electrons, and /Xc,^ is the matrix 
of spin velocities. Next, is an arbitrary function of 
the absolute value of electron momentum. This func- 
tion describes possible isotropic corrections to spin-orbit 
interaction, which may have the same origin as the non- 
parabolicity of the band spectrum. The case = I cor- 
responds to the p-linear spin-orbit Hamiltonian of the 
general form. The p-linear spin-orbit coupling terms 
appear in quantum wells due to both the structural in- 
version asymmetry^** (Rashba term) and bulk inversion 
asymmetry (Dresselhaus term) , the latter contribution is 
sensitive to orientation of the quantum well with respect 
to crystallographic axes. By solving the quantum kinetic 
equation for the matrix distribution function of electrons, 
with taking into account elastic scattering, the spin po- 
larization and spin currents are found on an equal footing 
and a relation between them is established. The calcu- 
lations not only provide analytical expressions for these 
quantities, but also demonstrate a need to reconsider the 
known results^'"^^ for static spin polarization based on the 
Hamiltonian (1) in special cases, when the symmetry of 
liafj allows zero spin precession for some chosen direc- 
tions of the spin vector. The examples of this kind are 
the quantum wells grown in [001] crystallographic direc- 
tion in the case of equal Rashba and Dresselhaus spin 
velocities and symmetric quantum wells (only the Dres- 
selhaus term is present) grown in [110] crystallographic 
direction. In particular, it is shown that the induced spin 
polarization in these systems remains zero even when the 
frequency of the applied field goes to zero. The calcula- 
tions also show that the static spin current in the electron 
systems described by the Hamiltonian (1) is zero even if a 
non-parabolicity of electron band spectrum is taken into 
account. The presented theory neglects the spin-orbit 
corrections to the scattering potential, which means that 
the effects of extrinsic spin-orbit coupling, such as the 
extrinsic spin currents and Elliot- Yafet spin relaxation, 
are not included in the calculations. 

The paper is organized as follows. In Sec. II we con- 
sider the quantum kinetic equation and present its an- 
alytical solutions. The expressions for the induced spin 
density vector and spin current tensor, obtained on the 
basis of these solutions, are given in Sec. III. In that 
section we also establish a general relation between these 
quantities and present a detailed analysis of two impor- 
tant cases, [001]-grown and [110]-grown quantum wells 
with both Rashba and Dresselhaus spin-orbit coupling. 
A relation between the induced spin density and electric 



current is derived and analyzed in Sec. IV. The obtained 
results and the limits of their applicability are discussed 
in Sec. V. 



II. GENERAL CONSIDERATION 



The Hamiltonian of the problem is written in the form 



H„ + hp + Vr 



where H„ is the Hamiltonian of free 



electrons in the absence of spin-orbit interaction, hp is 
the spin-orbit Hamiltionian, Vr is the potential of impuri- 
ties or other static inhomogeneities (the spin-orbit correc- 
tions to this potential are neglected, so only the intrinsic 
spin-orbit coupling is considered), and H^^* = — eEj-r is 
the Hamiltonian of the external perturbation due to the 
applied time-dependent electric field Ej (here e = — |e| 
is the electron charge). The calculations are based on 
the quantum kinetic equation for the Wigner distribution 
function ppt, which is a 2 x 2 matrix over the spin indices 
(see, for example, Refs. 36 and 37). For the spatially- 
homogeneous problem considered below, this equation is 
written in the form 



at 



h 



,ppt 



pt 



dp 



J{p\pt), (2) 



where J is the collision integral describing the elastic 
scattering. This integral is written below in the Marko- 
vian approximation and under the assumptions Uu) <^s 
and h/r <C e, where lu is the frequency of the applied 
perturbation, e is the mean kinetic energy of electrons 
and T is the characteristic scattering time. One has (see 
Ref. 37, problem 13.10) 



■^(PjP*) = 72 



dp' 

{2TTh)' 



rW^(lp-p'l) f 

J — OC 



(3) 



J(ep + hp)l' /h 



(A 



p't - ppt 



where W{\p\) = / dAre-^P'^"-/'^ ((Fr+ArK)) is the spa- 
tial Fourier transform of the correlation function of the 
scattering potential, Sp is the kinetic energy of electron 
in the absence of spin-orbit interaction (the energy spec- 
trum is isotropic but not necessarily parabolic), A +Q, 
and (p' ^ p) denotes the term obtained from the pre- 
ceding one in the square brackets by permutation of mo- 
menta. The integration over t' in Eq. (3) is carried out 
elementary, but the resulting expression is rather lengthy 
and, for this reason, is not presented here. 

Searching for the linear response to the Fourier com- 
ponent Ee~"^* of the applied electric field, we repre- 
sent the matrix distribution function in the form ppt = 
^^q) _|_ ^^g-iwt^ where fp is the Fourier component of 
the non-equilibrium part of the distribution function and 



(eq) 



is the equilibrium distribution function, 

f^l) — — [/£p_|_^Qp -|- fep—HQp] 



+ - 



fir 



20. 



'[fep+HQp /ep-ROp]) 



(4) 
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which is expressed through the Fermi distribution f^. 
Here and below, Hp = \ flp\. The function (4) does not 
lead to spin polarization of electron system because its 
matrix part is antisymmetric in momentum. Since a sub- 
stitution of the function (4) into Eq. (2) makes both the 
commutator and the collision integral (3) equal to zero, 
one has a closed integral equation for /pi 



«w/p + - hp, fp 



+ eE- 



dp 



= Jif\p), (5) 



which determines the linear response of the electron sys- 
tem. 

The induced spin density is defined as 
I f _dp 



s = 



jTr{&fp), 



(6) 



and the induced (non-equilibrium) spin current density 
is given by the tensor 



rTr({(Ta,%(p)}/p), 



(7) 



where u(p) = d{ep+hp)/dp is the group- velocity matrix, 
Tr denotes the matrix trace, and {a,b} = [ah + ba)/2 
denotes the symmetrized matrix product. The expression 
(7) describes the flow of the spin polarized along a in the 
direction 7. One may also introduce the average spin S = 
s/«2D, where n2D = (27r?i)~^/ dpTrpp is the electron 
density. The tensors of spin polarizability, Xa/3) and spin 
conductivity, S"^, are introduced according to 

Sa = Xa0{uj)E0, = T.^0{uj)E0. (8) 

It is assumed in the following that the spin-splitting 
energy 2hVlp is small in comparison to the mean energy 
of electrons. Then it is convenient to apply an efficient 
method of solution of Eq. (5) based on the expansion of 
the collision integral in series with respect to the small 
parameter UClp/e; see Refs. 5, 25, 31, and problem 13.11 
in the book 37. Using the spin-vector representation 



fp 



& ■ fp and retaining only the terms of the 



first order in Jlp under the collision integral, we obtain 
coupled equations for scalar and vector parts of the dis- 
tribution function: 



1 TP ^[/gp-t-ftOp + fe 

2^ ' dp 



27r 



^^H^(|p-p'|)[(fO,-fXv-e,) 



-h{ilp - Up,) . (fp, - fp) 



dd{ep' - £p) 



de. 



p' 



(9) 



and 



- iuiip + eE 



27r 



dp 
dp' 



2Vtr 



M/(|p-p'|)[(fp,-fp) 



2["p X fpj - ^ J 



p' 



.(10) 



The expressions containing formal derivatives of the 5- 
functions under the integrals should be evaluated using 
integration by parts. The Fermi distribution functions 
standing in the field terms also can be expanded in series 
of Op . Then one can see that the iterational expansions 
of fp and fp start with the terms of zero and first order 
in Clp, respectively. For this reason, the last term under 
the collision integral in Eq. (9) can be neglected, and the 
solution of this equation is 



eE-Vp df,^ 



,(1) 



iw ^^p 



(11) 



where Vp = de-p/dp is the group velocity of electron in 
the absence of spin-orbit interaction. Here and below, the 
relaxation rates appearing in the problem are defined as 



2tt 



dp' 

{2Trny 



■W{\p-p'\)[l-cos{ne)]6{ep,-ep), 

(12) 

where 9 denotes the scattering angle pp'. The rate Up"'^ 
describes relaxation of the n-th angular harmonic of the 
distribution function. 

Equation (11) describes the Drude response of the elec- 
tron system. The next correction to fp is of the order 
of {hflp/e)"^. This correction is essential for calculation 
of the frequency-dependent conductivity and dielectric 
function of 2D electrons with spin-orbit splitting'^*''^^ (see 
Sec. IV), but it is not important for calculation of the 
induced spin polarization and spin current. After sub- 
stituting the expression (11) into the last term of the 
collision integral in Eq. (10), this term is unified with 
the field term on the left-hand side of Eq. (10). As a 
result, one gets a closed equation for the vector- function 
fp: 



-iu}fp-2[np X fp]-hFp(w) 



= ¥ /(^^(|P-P'l)(fp'-W^f'-^^')- (13) 

The vector F contains both isotropic and anisotropic con- 
tributions: 



Fp(a;) = eflE^i^H'^) " e(E • p)np^W(w), 



(14) 



^«(a;) = Rpiu;){l^p - tu) -p^^M^.i^), (15) 



^W(a;) 



2h 



dp^ 

^.(*-P^^(-f--) 



,(16) 



where JIe is the constant vector with components 



^p'' - iuj 9ep 



(17) 



and 



-^P ^P 2 dep dp^ J ^ 
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is a relaxation rate. Note that Vp goes to zero in the Umit 
of short-range scattering potential. 

It is convenient to search for the solution of Eq. (13) 
in the form 



e(E.p)np:rW(^) 



(2) 



gp- 



(19) 



Since the first term of this expression is proportional to 
rjp, it does not contribute to the vector product in Eq. 
(13), thereby representing a non-precessing part of the 
sohition. Note that the angular average of this term is 
directed along 17e. The substitution (19) leads to the 
following equation for gp: 

- iujgp - 2[rjp X gp] + iujGp{uj) 

= ^ ^m\p - p'i)(gp' - gp)]|p'i=iPh (20) 

where rup = | {dp^ / dsp) is the p-dependent effective mass 
and 



Gpito) 



rip 



e^p{v^p 



(2) 



{Vp^ — iij) {Up^' — iij) 



(2) 



(21) 



The collision integral in Eq. (20) is already reduced, 
by means of integration over the absolute value of p', 
to the integral over the angle (p' of the vector p'. Ow- 
ing to the substitution (19), the inhomogcncous (field- 
dependent) term of Eq. (20), iu)Gp{u)), is proportional 
to the frequency ui and isotropic in the momentum space. 
With the aid of the definition (21), it is convenient to 
write the angular-averaged distribution function fp = 



(27r) ^ Jq'' dipfp as 



-Gp{oj) + Qp{u)+g 



(22) 



Qp{u>) = -e^E 



d\fRp{u 



We also point out the exact relation 



A. Short-range scattering potential 

Let us consider the limit of short-range scattering po- 
tential, when W{\tp — p'|) is replaced by a constant W 

and t'p"^ = Vp = rUpW/h^ for any number n. The mo- 
mentum dependence of the scattering rate is associated 
with possible non-parabolicity of the band spectrum and 
has to be ignored in the parabolic approximation. Since 
the right-hand side of Eq. (20) is reduced in this case to 
Vp{gp — gp), a regular way of solving exists. The solution 
is 



gp 



{{up - iw)'^Ap 



Vpi^p - iw)A2 d£p 
+2{vp - ia;)[np X Ap] + 4J2p(J2p • Ap)} , (24) 

where Ap = {vp — iiv)^ + 4f2p and Ap is the vector with 
components 



= eUTp-'U^^. 



(25) 



Here Tp ^ denotes the matrix inverse of the symmetric 
3x3 matrix 



= 4 (f]«f2g - do^0ni)/Al + S^- (26) 



obtained as a result of angular averaging. The whole 
solution, according to Eqs. (16)-(19), is 



f, = -e(E.p)n,||, 



Up — iij dsp 



+ gp, (27) 



and its angular average is written as 

1 



Vp dsp 



lUJ 



e^pflE Ap -|-Qp(a;). 



(28) 



This vector determines the magnitude and the direction 
of the induced spin polarization. 



fp = - ^[fip X gp] -I- Qp(a;), 



(23) 



which is obtained by applying the procedure of angular 
averaging to Eq. (20) and by using Eq. (22). 

In spite of the isotropy of the term (21), Eq. (20) 
requires a numerical solution. The physical reason for 
this is the effect of precession in the presence of angular- 
dependent scattering. The angular dependence of the 
vector product [Jlp x gp] , in the general case, is different 
from that of gp standing there, and the standard method 
of solution, based on expansion of the distribution func- 
tion in series of angular harmonics, leads to an infinite set 
of coupled equations. There are, however, a number of 
important situations when Eq. (20) can be solved analyt- 
ically. These situations are described in the subsections 
below. 



B. Isotropic spin splitting 

The next exactly solvable situation is realized when 
the energy spectrum of electrons remains isotropic in the 
presence of spin-orbit coupling. In other words, fJp — fip 
depends only on the absolute value of p. This imposes 
certain constraints on the matrix 



l^xx^xy ^ f-^'yxf-^'yy ^" l^zxl^zy 0- 



(29) 



This situation is realized, for example, in [001]-grown 
quantum wells {x ||[100], y || [010], z || [001]) with only 

Rashba or only Dressclhaus type of spin-orbit coupling 
and in [lll]-grown quantum wells [x \\ [112], y \\ [110], z \\ 
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[111]) with both types of couphng. The solution is writ- 
ten in the form similar to that of Eq. (24). By introduc- 
ing the isotropic quantity = {up^^ —iui) {v^^ —iui)+4:Qp 
and the vector 



(30) 



where Up ^ denotes the matrix inverse of 



uf = 4{n»n^p-s^pnl)/Al + 5^0^, (3i) 



we find 



('2") 

iuj{iyp ' +Vp- iuj) dfe^ 



gp = 



x{(t.W-iu;)(./f -iu;)Bp (32) 

+2(i/f ) - iuj)[np X Bp] + Aftpiftp ■ Bp)} . 

It is easy to see that the results (24) and (32) become 
equivalent if one assumes isotropic spin splitting in Eq. 
(24) and short-range scattering in Eq. (32). The angular 
average of the whole solution is 

= [2), (1) I «?p"e - —Bp I +Qp(a>). 



(33) 



C. Fixed precession axis 



There is also a special case, when Eq. (20) is solved 
in the most simple way. This happens when the vector 
product [rip X He] is zero for arbitrary p and E. In other 
words, the symmetry of the matrix jiap should allow ex- 
istence of a fixed (momentum-independent) precession 
axis. This imposes the following constraints: 



(34) 



Under these conditions, both Sip and rie arc directed 
along the fixed precession axis, without regard to direc- 
tions of p and E. The examples are [001]-grown quantum 
wells {x II [100], y \\ [010], z \\ [001]) with equal absolute 
values of Rashba and Dresselhaus velocities, when the 
precession axis is in the quantum well plane at the angle 
of 7r/4 or — 7r/4 with respect to the main crystallographic 
axes, and [110]-grown wells {x || [TlO], y \\ [001], z \\ [110]) 
with only Dresselhaus type of coupling, when the preces- 
sion axis is perpendicular to the quantum well plane. The 
solution of Eq. (20) in this case is non-precessing (makes 
the vector product equal to zero) and isotropic (makes 
the collision integral equal to zero): 



(35) 



This vector is directed along He- The averaged whole 
solution fp is also directed along Sl^- 



Qp(w). 



(36) 



D. Static limit 



If the frequency u) goes to zero, the solution of Eq. 
(20) is trivial, gp = 0. Therefore, the fimction fp from 
Eq. (19) with gp = and a; = describes the static 
spin-dependent response in the general case. The only 
exception is the special case considered in the previous 
subsection, when there exists a non-zero solution, gp = 
Gp(0). The averaged distribution function for the static 
limit is 

for the general case. In the special case only the last term 
of this expression remains. 



III. SPIN RESPONSE 

To describe the spin response, one should calculate the 
integrals over momentum p in Eqs. (6) and (7). It is 

convenient to separate the angular averaging from the 
integration over the squared absolute value of momen- 
tum, p^, according to rfp = | dp^ dip. Then, after using 

the representation /p = fp + o" • fp and taking the ma- 
trix trace, the density of the induced spin polarization is 
given by 



f 

Jo 



dp^ 



_f 



(38) 



and the density of non-equilibrium spin current is 



f 

Jo 



dp^ 



+Ma/3fp (^/37^P + 2(a^p/ap2)p^p^) ]. 



(39) 



The second term, which appears in the expression (39) 
owing to the spin-orbit correction to the group velocity, 
gives zero contribution because, according to Eq. (11), 
the scalar part of the distribution function is antisymmet- 
ric in momentum. If the frequency lu is zero, the vector 
part of the distribution function is symmetric in momen- 
tum, so the first term of the expression (39) also gives 
zero contribution. Therefore, the non-equilibrium static 
spin currents do not exist for the model described by the 
spin-orbit Hamiltonian (1). At non-zero frequency, the 
spin currents are associated with g-contribution to the 
distribution function, because the first term of Eq. (19) 
is symmetric in p: 



dp"^ 



-gpP7> 



(40) 



where = {qj,q^,qj)- Looking at the expressions (24) 
and (32), one can conclude that it is the second terms in 
the braces of these expressions that are responsible for 
the spin currents. 
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On the other hand, the induced spin density is deter- 
mined by the angular-averaged symmetric part of fp and 
exists in the static regime as weh. Applying either the 
exact relations (22) and (23) or the expressions (28), (33), 
(36), and (37) describing different physical situations, one 
should always ignore the term Qp{uj) which represents a 
full derivative over p"^ and, for this reason, does not con- 
tribute to the integral (38). The static spin polarization 
obtained in this way is 



47r/i ./n ^ dSn 



1 

72)" 



241)42) de. 



(2) ^^P 



(41) 



In the parabolic approximation Sp = /2m and at = 1 
this expression is rewritten as 



Hi 



27r?i 



where ri"'* = l/vp''' are the relaxation times. For degen- 
erate electron gas, when {dfs/de) ~ —5{e — ep) and ep 
is the Fermi energy, the integral over energy is taken in a 
straightforward way. In the case of short-range scattering 



de 



de 



^(2) , ^(Di^lnr, 



(2) 



2 dine 



(42) 



point of view, this paradox is related to the fact that 
in the static limit the additional part of the distribution 
function, gp = Gp(0), which cancels the contribution of 
the first term in the expression (19), still exists, only in 
the case of a fixed precession axis. One can say that 
the dependence of the induced spin polarization on the 
parameters (components /^a/s) of the spin-orbit Hamilto- 
nian is non-analytic in the region of parameters where 
the fixed precession axis appears. This means that the 
result for the spin polarization depends on the order of 
limiting transitions. If first co is aimed to zero and then 
the precession axis is fixed, the polarization is finite. If 
first the precession axis is fixed and then uj is aimed to 
zero, the polarization is zero. More details on the issue 
of non-analyticity will be given below in this section, by 
considering concrete examples. It is important to state 
that the consideration given above docs not provide spin 
relaxation mechanisms other than the DPK mechanism. 
Inclusion of the Elliot- Yafet mechanism (see Refs. 42 and 
43 for the 2D case) can lead to a finite relaxation of the 
spins oriented along the fixed precession axis and, there- 
fore, to a finite induced spin polarization for this special 
case; see Sec. V for more discussion. 

Using Eq. (23) together with Eqs. (20), (38), and 
(40), one can obtain an exact relation connecting the 
induced spin polarization and spin current. From Eqs. 



potential (ri") = r = 1/i.) one arrives at the well-known g^) and (23) one has s = -{27riu) ^ dp^i^^ x gp]. 

Smce [SZp x gpj = x gpJp/3, where /x^ is the vector 

with components /i„/3) we find, for the case of parabolic 
band and ^p = 1, 



result^^ Sa = XajsEfs with Xap = -fiap[emT/2nh'^] 
valid for any linear spin-orbit Hamiltonian including the 
case of Rashba spin-orbit coupling studied in the early 



papers 
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It is important that Eqs. (41) and (42) are not valid 
in the special situation when a fixed precession axis ex- 
ists; see subsection C of Sec. II. In fact, a straightforward 
substitution of Eq. (36) into Eq. (38) shows that, at arbi- 
trary frequency of the applied field, the spin polarization 
does not appear in this special situation. The physical 
explanation of this remarkable property is based on the 
facts that in the case of a fixed precession axis (a) there is 
no spin relaxation^" by the Dyakonov-Perel-Kachorovskii 
(DPK) mechanism^i for the spins directed along this axis 
and (b) the induced spin polarization can, in principle, 
appear only along this axis, without regard to the direc- 
tion of the applied electric field. Imagine that the electric 
field is abruptly turned on. The anisotropic distribution 
of electrons over momenta, which determines the electric 
current, is established during the momentum relaxation 
time l/i>p^\ However, the distribution of electrons over 
spins, which determines the spin density, is established 
during the spin relaxation time, and the corresponding 
transient process becomes infinitely slow if the spin re- 
laxation is absent. Therefore, if a periodic alternating 
field acts on electrons in a sample with a fixed preces- 
sion axis, the spin density cannot react to the field at 
any frequency u, and the spin polarizability is zero. The 
absence of the static spin polarization in the case of a 
fixed precession axis cannot be revealed by considera- 
tion of the static spin response alone. From the formal 



2m 
iJiuj 



(43) 



This relation is remarkable in the sense that it does not 

contain any parameters describing scattering. The re- 
lation (43) is also valid for a non-parabolic model pro- 
vided that the momentum dependence of l/nip and ^p 
is the same; in this case the effective mass m should be 
replaced by the momentum- independent quantity rrip^p. 
Using Eq. (43), one can directly relate the low- frequency 
behavior of the spin current to the static induced spin 
polarization. Equation (43) can be derived from the 
continuity equation (for the case of Rashba model see 
Ref. 44 and Sec. 64 of Ref. 37). In detail, when 
Eq. (2) is summed over the momentum p, the collision- 
integral contribution and the electric-field term vanish, 
and Eq. (43) directly follows from the relation between 
the momentum-averaged first and second terms on the 
left-hand side of Eq. (2). This observation also demon- 
strates that the validity of Eq. (43) is not restricted 
by the assumption of linear response used in this paper: 
the non-equilibrium part of the distribution function p 
is not necessarily small. On the other hand, Eq. (43) is 
specific for the chosen form (1) of the spin-orbit Hamil- 
tonian. Applying Eq. (43) to the model including both 
Rashba and Dresselhaus linear spin-orbit coupling terms 
in [001]-grown quantum wells (see subsection A below), 
one obtains the relations between the components of the 
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spin density and spin current recently derived"'^^'^^ for this 
model. The corresponding rclations'^''''*^ for the Rashba 
model also follow from Eq. (43). 

Since it is established that the non-parabolic correc- 
tions to the band spectrum and the corrections to the 
spin-orbit Hamiltonian do not lead to qualitative modifi- 
cations of the induced spin polarization and spin current, 
we neglect these corrections in the following, by assum- 
ing Ep = /2ra and = 1. Let us consider first the fre- 
quency behavior of the induced spin polarization and spin 
current in the quantum wells described by the Rashba 
model. The non-zero components of the spin-velocity 
matrix are fixy = — Mys = ^-Rj where is the Rashba 
velocity. The spin splitting described by the Hamilto- 
nian (1) is isotropic in this situation. Using Eq. (33) and 
taking into account that the matrix (31) is diagonal, we 
obtain 



em _ , , 



(44) 



where the frequency-dependent function 




Sp dvf^ jdsp 
9 (1) 

^ Up ' ~ ILO ^ 



2{vp/nf 



2[vvlWvf^- i^\{vv'~ ioj){i^p'~ *w) + 4.{vp/n) 

has dimensionality of time. Equations (44) and (45) gen- 
eralize the result of Ref. 31 obtained in the limit of short- 
range scattering potential to the case of arbitrary scatter- 
ing potential. They also describe [lll]-grown quantum 
wells with both Rashba and Dresselhaus spin-orbit cou- 
pling, where fi^y = —f^yx = vr + 2v£)/V3 and vd is the 
Dresselhaus velocity. In this case, vn in Eq. (44) should 
be merely replaced by vji + 2vd/V3 and the direction of 
the spin polarization vector s with respect to E is the 
same as for the pure Rashba coupling. In particular, 
s = [emvT(LU,v)/2Trh'^](—Ey, Ex,0) where v is cither vj^ 
or vn + 2vd/V^, respectively. Next, in symmetric [001]- 
grown quantum wells, where only the Dresselhaus spin- 
orbit coupling is present and the non-zero components 
are fj,yy = —jJixx = vd, Eq. (44) with T({jj,vd) is also 
valid, leading to s = [emvDT{(jj,VD)/ 21^11 ]{Ex,—Ey,Q). 
In all these cases, the spin polarization is in the quantum 
well plane [sz = 0), and the direction of this polariza- 
tion is frequency-independent. The spin currents can be 
either calculated directly or extracted from Eq. (43). 
There exist only the currents of ,2-polarized spins, and 
these currents are directed perpendicular to the applied 
field: 



(1) 



(2) 



-(45) 



Tx 
1y 



ieujT{u;, v) 



Ex 



(46) 



where T{lo,v) is given by Eq. (45) with v = vr for the 
pure Rashba coupling and v = vii-\- 2uo/\/3 for [111]- 

grown quantum wells. For symmetric [001]-grown quan- 
tum wells one should use v = vd and change the sign of 



the right-hand side. Therefore, the symmetry properties 
and frequency dependence of the spin currents remain 
the same for all important cases of isotropic spin split- 
ting considered in this paragraph. In the limit of short- 
range scattering potential and in the case of degenerate 
electrons, Eqs. (46) and (45) with v = vr give the result 
obtained previously in Ref. 16. The universal behavior ^'^ 
of the spin-Hall conductivity E^^ : 



Yjy^ exists in the 



limit vpF u) '> UpJ , where pp is the Fermi momen- 
tum. 

Now we turn to more complicated situations when the 
anisotropy of spin splitting is essential due to combined 
effect of both Rashba and Dresselhaus spin-orbit cou- 
pling. These cases are considered in the following sub- 
sections in the limit of short-range scattering potential, 
when the expressions obtained in subsection A of Sec. II 
are valid. 



A. [001]-grown quantum wells 

Let us study the case of [001]-grown quantum wells. 
If the Cartesian coordinate axes are chosen along the 
principal crystallographic directions, there are four com- 
ponents of the spin- velocity tensor: 



f^xy 



l^yx — ^i?, ^yy 



VD- 



(47) 



The spin splitting depends on the angle (f of the vector p 
according to Op = {v'^+vfj){p/?i)^-2vRVD{p/?i)'^ sm{2ip) 
and Ap = a — bsm{2ip), where 

a={u- iujf + A{v% + vl)ip/nf, b = 8vRVD{p/hf. 

(48) 

The angular averaging in Eq. (26) results in ^ = T^^ = 
{ly - iujf/2r - 1/2 + iuj/v and T^^ = T^f = b/2r - 

(ajr — l)(w|. + vj-i) / Av rv D ^ where r = \J — b-^. After 
some transformations, we obtain the relation defining the 
tensors of spin polarizability and spin conductivity: 



em 



2'!rh- 



-I 

V Jo 



Kp{LO) 
Kp{u)) 



Dp{LO) 

-Kp{Lu) 
-Kp{u)) 



dsp 

Ex 
E,. 



(49) 



and 



IBLO V% 



A-Khv Avrvd Jo Dp{uj 

9p{uS) -9p{oj) 
6p{uj) -6p{u)) 

where the denominator is given by 



dfe, 

dsp 

Ex 
Ell 



(50) 



Dp{^) 



(- 

Vr 



- 1 



V 



1 + - 



(51) 
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The elements of the matrix in Eq. (49) are 



4vd 



24 



V 



(52) 



and Kp(a;) is obtained from this expression by the per- 
mutations Vfi vd- The elements of the matrix in Eq. 
(50) are given by 



W = 1 



V 



and 



V r 



(53) 



(54) 



Equations (49) and (50) demonstrate the symmetry re- 



= -E: 



!/!/' 



and 



lations Xxx — Xyy^ Xxy — Xyxj 

yz — —yz 

xy yx' 

Though the matrix in Eq. (49) retains the symme- 
try of jJiaffi the ratio Kp{uj) / Kp{u)) is not equal to vr/vd- 
For this reason, the direction of spin polarization in the 
plane is different from the direction of Oe and depends 
on the frequency. The direction of the spin current is 
not perpendicular to the direction of the field and is also 
frequency-dependent. The components of the spin den- 
sity and spin current are related according to Eq. (43), 
which can be written, for this particular case, in the 
fQj.j^i8,32 _ (^ifii_^/2m){vRSx + VDSy)/{v\ — vjj) and 

= (ihuj/2m){vDSx + VRSy)/{v\ - v^,). 

In the collisionless limit, ^ 0, Eqs. (49) and (50) are 
reduced to the results obtained in Ref. 32. Note that the 
formal substitution z/ — > makes the spin currents finite 
at ti) ^ because the denominator Dp{uj) is reduced 
to —{lo/v)'^, while the functions (53) and (54) become 
proportional to u/v. The static spin currents for the 
systems described by the spin-velocity matrix (47) have 
been also studied in the collisionless limit in Refs. 45 
and 46. All these studies show that, as the ratio vr/vd 
is varied, the spin current reverses its sign going through 
zero at v\ = vj^. This general property, also reflected in 
Eq. (50), follows from the Berry phase analysis.^''' 



In the case ^ 



as already mentioned, there exists 



a fixed precession axis directed at the angle of 7r/4 (or 

— 7r/4) in the quantum well plane. Therefore, not only 
the spin current, but also the induced spin density goes 
to zero in this case. This behavior is demonstrated in 
Figs. 1-3, where the calculated real parts of the compo- 
nents of spin polarizability and spin conductivity tensors, 
ReXa/9 and ReS^^, are plotted as functions of the ratio 
of Rashba and Dresselhaus velocities. The case of de- 
generate electron gas is assumed. The spin polarizability 
is expressed in the units of static polarizability for sym- 
metric [001]-grown quantum wells, xo = emvD/'^T^Ti'v- If 
the frequency decreases, both Xxx and Xyx, as expected, 
approach their static values xo and {vR/v]j)xa: respec- 
tively. However, this never happens at vr = vd, when 
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FIG. 1: Dependence of the real part of spin polarizabil- 
ity in [001]-grown quantum wells on the ratio vr/vd at 
vdPf/Tii' = 1 for several frequencies. The polarizability is ex- 
pressed in the units of xo = emvo/'i'K'h'^v. The diagonal and 
off-diagonal components are shown by the solid and dashed 
lines, respectively. 




FIG. 2: Dependence of the real part of spin condxictivity in 
[001]-grown quantum wells on the ratio vr/vd at vdPf /hv = 
1 for several frequencies. The spin conductivity is expressed 
in the universal units. The diagonal and off-diagonal (Hall) 
components are shown by the solid and dashed lines, respec- 
tively. The difference between them is negligible at small oj. 



the polarizability remains exactly zero at arbitrary uj. 
The real part of the spin conductivity (Fig. 2) shows 
prominent peaks near the point vr — vd at small u), 
though its behavior is analytic. The imaginary parts 
of Xxx and Xyx (not shown) also have sharp peaks in 
the vicinity of vr = vd at small cj. The depression 
of the spin polarizability in the region vr ~ Vd is ex- 
tended with the increase of the disorder, as shown in 
Fig. 3. This also means that, especially for the "dirty" 
case vdPf hv, the frequency dispersion of the spin po- 
larizability remains significant at very small frequencies 
if Vr is close enough to vjj. The corresponding behavior 
is illustrated in Fig. 4 and is explained by the reduc- 
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tion and disappearance of the DPK spin relaxation as vr 
approaches to vd. 




FIG. 3: Effect of disorder on the real part of spin polar- 

izability. The "clean" and the "dirty" cases correspond to 
VDPF/hu = 10 and voPF/hv = 0.1, respectively. 



Applying the equations listed in subsection A of Sec. 
II, we find that the field can induce both y- and 
^-polarized spins, while the field Ey induces only x- 
polarized spins: 



emvuEy 
s-r = -TT-^ / de 



2TTh^v Jo ""^ V 



K-{uj) 



emvRE. 



and Sz = —{vD/^VftjSy, where 




and 



(57) 



c — d \ iu> 

— -7-1 +— , 59 
c + d / V 
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FIG. 4: Low-frequency dispersion of the real part of spin 
polarizability in [001]-grown quantum wells in the vicinity of 
the fixed precession axis {vr is close to vd)- 



B. [110]-grown quantum wells 

The next case wc consider is [110]-grown quantum 
wells. If the Cartesian coordinate axes are chosen with 
OZ perpendicular to the quantum well plane and OY 
along the principal crystallographic direction {x \\ [TlO], 
y II [001], z II [110]), there are three components of the 
spin- velocity tensor: 



Mx!/ = — Mycc = Vr, IJ.ZX = 



(55) 



One has fip 



(4 + vl/8){p/hr + K/8){p/hr cos(2^) 
and Ap = c + dcos{2(p), where 

c={u- iujf + {Avl + vl/2){p/hf, d = vn{pinf/2. 

(56) 



The component is a symmetric function of vr, while 
Sx and Sy are antisymmetric functions of vr. 

The spin currents appear for y- and z-polarized spins 
and flow in the direction perpendicular to the applied 
field: 



qy = 



,,2 /„,2 7T\^x, Qy — 



rnvDi^vl/vl + 1) 



ihujVD 



(61) 



and 



q', = {2vR/vD)ql q'y = {2vR/vD)ql (62) 



The relations (61) and (62) satisfy the general require- 
ment (43). The components q^. and qy arc symmetric 
functions oivR, while q^. and are antisymmetric func- 
tions of Vr. 

For symmetric quantum wells, where vr = 0, both the 
polarization and the spin currents disappear for arbitrary 
LO. If both Vr and w go to zero, the ^;-component of the 
spin polarizability tensor, Xzx{^), is described by a sim- 
ple formula applied to the case of degenerate electrons: 



ZX — 



Xo {2vr/vdY^f 
2 {2vr/vdYif - iw' 



(63) 



where — \/v'^ + (vdPf /K)"^ — v. Equation (63) illus- 
trates the non-analytic behavior of the spin polarizability 
in the vicinity of the fixed precession axis. The calculated 
dependence of the real part of Xzx on the ratio vr/vd, 
shown in Fig. 5 for the case of degenerate electron gas, 
is similar to the dependence of Xxx shown in Fig. 1 for 
[001]-grown wells. The region of depression near vr = Q 
is extended in the "dirty" limit v^pp <^ hv, as it is seen 
directly from Eq. (63). The behavior of the component 
Xyx remains analytic at vr ^ and w — > 0, though it 
is strongly affected in the vicinity of vr = at low fre- 
quencies. In contrast, the component Xxy (not shown in 
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Fig. 5) stays close to its static value — (f_R/fz5)xo up to 
the frequency region lu ^ v. The real part of the spin 
conductivity shows peaks in the vicinity of vr = 
at low frequencies, while the low-frequency behavior of 
J^^y is monotonic in this region, see Fig. 6. According 
to Eq. (63), the frequency dispersion of the spin po- 
larizability remains significant at very low frequencies if 
vr is small enough. It is remarkable that in the "dirty" 
limit the characteristic rate describing this dispersion, 
{^vn/vof^p, is equal to 2{viiPF/K)^ /f, which is the 
DPK spin relaxation rate for the Rashba model. 
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FIG. 5: Dependence of the real part of spin polarizability in 
[110]-grown quantum wells on the ratio vr/vd at vopp/fii' = 
1 for several frequencies. The polarizability is expressed in 
the units of xo = ernvnl'^T^'f^v. 
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FIG. 6: Dependence of the real part of spin conductivity in 
[110]-grown quantum wells on the ratio vr/vd at vopp/'hi' = 
1 for several frequencies. 



IV. CURRENT RESPONSE 



to the frequency dispersion of the conductivity (or di- 
electric function) of 2D electron layers with spin-orbit 
interaction. Some relations of this kind have been estab- 
lished previously"^^'"^^ in the collisionless approximation. 
In this section the corresponding relations are obtained 
and analyzed for the general p-linear model of spin-orbit 
coupling described by the Hamiltonian (1), with taking 
into account the electron-impurity interaction. The cal- 
culation is based on the kinetic equation (2). The electric 
current density is defined as 



dp 



rTr(u(p)/p). 



(64) 



Below we neglect the non-parabolicity of the energy spec- 
trum and the deviation of the Hamiltonian (1) from the 
linearity, by putting nip = m and = 1. Let us multi- 
ply Eq. (2) by ep/m, sum it over p, and take the matrix 
trace. Using Eqs. (6) and (64), one has the exact relation 



e^Ef3n2D 



m 



f d] 
= V(^ 



dp p/3 



,^Tr^(/|p). 



(65) 

The right-hand side of this equation should be set at 
zero in the collisionless approximation. To consider the 
collision-induced contribution in the general case, one 
should calculate the distribution fmiction fp with the 
accuracy up to the terms ^ {hilp/s)'^. Instead of doing 
this, we consider the limit of short-range scattering po- 
tential, when the right-hand side of Eq. (65) is exactly 
transformed to — z^j/j. Therefore, 



J0 



m(z/ — iuj) 



2iuj 



(66) 



Equation (66) contains the usual Drude term and the 
term induced by the spin-orbit interaction. Together 

with Eq. (43), this equation establishes the relationship 
between the electric current, spin current, and induced 
spin density. The validity of both Eq. (43) and Eq. (66) 
is not restricted by the assumption of linear response. 

Expressing the currents through the electric conductiv- 
ity CFi3x{ijj) and spin conductivity S^^(t<;), one can write 
the equation that relates these tensors: 



. e n2D , (SO) 



Sso) 



m{v — iuj) 
4me 



n{u - iuj) 



(67) 



In general, the spin-orbit term a^^^^ brings non-diagonal 
contributions to a^x. This should lead to a weak Hall 
effect in the absence of magnetic field at finite frequen- 
cies. For the Rashba model, when only the components 
IJ'xy = —f^yx = vn exist, the conductivity is diagonal and 



isotropic, ' a 



47 JSO) 



/3A 



(^^^^^Si3x, where 



It is important to relate the behavior of spin polar- 
ization and spin currents studied in the previous section 



(SO) ^ 4me4 



(68) 
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In the coUisionless hmit, this equation gives a relation^'' 
between the imaginary part of cr*-'^'^-' and the real part of 
the spin-Hall conductivity S^^^. 

In the case of [001]-grown quantum wells with both 
Rashba and Dresselhaus types of coupling, 



(SO) 



a 



xX 



fi{v — iuj) 



(69) 



The non-diagonal part of o"^^'^'' is related to the diagonal 
components of the spin conductivity, while the diago- 
nal part is expressed through the spin-Hall conductivity 

(such an expression has been recently established^"^ in 

the coUisionless regime). The tensor a^^^'' can be diago- 
nalized by in-plane rotation of the Cartesian coordinate 
axes OX and OY . In contrast to the spin polarization, 

the quantity (J^p^^ is analytic at v\ = vljj. 

In the case of [110]-grown quantum wells the conduc- 
tivity is diagonal (in the chosen coordinate system) but 
anisotropic: 



"yy 



Ameiv 



--D 



/4) 



Ti{v — iui) 



yx 



(70) 



The expressions for the components S^^ entering Eqs. 
(69) and (70) are obtained from the expressions for spin 
currents presented in subsections A and B of the previous 

section. 

Although the contribution cr'f^'' is small as {hflp/e)'^ 
with respect to the Drude conductivity, its frequency de- 
pendence has qualitatively new features. If fip ^ ly, 
so that the spin-split states are well-defined, the term 
describes resonance absorption of electromagnetic 
radiation, typically in the THz region, associated with 
transitions between these states. In the case of isotropic 
spin splitting^^ (the Rashba model is considered below) 
and degenerate electron gas, the resonance takes place^^ 
a.t LO uir = 2\vn\pF /Ti- Substituting in Eq. (68) the de- 
tailed expression for the spin conductivity, see Eq. (46) 

and Eq. (45) with Up^^ = v^^ = u, one has an equation 



(71) 



which describes the resonance imder consideration at 
ujr ^ i^- Of course, this resonance also exists in 
the frequency dependence of the spin conductivity and 
spin polarizability.'^^ Another important feature following 
from Eq. (71) is the presence of low-frequency dispersion 
under the opposite condition, w^. <C v- This dispersion 
appears when w is comparable to the DPK spin relax- 
ation rate. Since this rate, for in-plane spin polarization, 
is given by vdp — 1{vji-pF jfif' jv, Eq. (71) in the limits 
<C I' and w <C gives 



(72) 



'DP 



Therefore, cr^'^'^' essentially depends on lj in the region 
of frequencies w <C J^, when the Drude conductivity still 
remains frequency- independent. Another contribution 
to the frequency dependence of the conductivity in this 
region exists owing to weak localization.'*® Though the 
weak-localization correction is larger in magnitude than 
RecT^'^'^^ its frequency dependence is slow (logarithmic), 
and can be distinguished from the dependence given by 
Eq. (72). 



V. SUMMARY AND DISCUSSION 

In this paper, the electric-field-induced spin density 
s and intrinsic spin current q/3 in 2D electron layers 
described by the general p-linear spin-orbit interaction 
Hamiltonian (JaUa/sPp = ^ ■ fJ-pPff are studied in the 
classical region of frequencies u>. The consideration is 
done for macroscopic systems and at zero magnetic field. 
The quantities s and q/3 are closely related to each other 
through Eq. (43) following from the balance equation 
for spin density. To find them, a careful analysis of the 
quantum kinetic equation is carried out taking into ac- 
count interaction of electrons with impurities or other 
static inhomogeneities. The presented results are valid 
for arbitrary correlation between the spin splitting en- 
ergy 2?iOp, disorder-induced broadening hu, and energy 
hu), under condition that all these energies are small 
in comparison with the characteristic kinetic energy of 
electrons. A complete analytical solution of the linear- 
response problem for the matrix distribution function is 
given in the limit of short-range scattering potential (see 
subsection A of Sec. II). This solution is applied to [001]- 
and [110]-grown quantum wells with both Rashba and 
Dresselhaus types of spin-orbit coupling (subsections A 
and B of Sec. HI). All other situations when analyt- 
ical solutions exist are described is subsections B, C, 
and D of Sec. II. The theory also takes into account 
the isotropic energy-dependent corrections to the effec- 
tive mass rn (non-parabolicity effect) and to the spin- 
velocity matrix fiap ■ This is reflected by the substitutions 
m — > nip and ^ai3 — * ^plJ^ajS assumed from the beginning 
of the consideration. Since these weak corrections do not 
lead to qualitative effects (in particular, it is shown that 
the static spin currents remain equal to zero in the pres- 
ence of these corrections), they are ignored in the most 
part of the applications, starting from Eq. (42). 

The main approximation of the present consideration 
is the neglect of the spin-dependent contribution to the 
scattering potential. This contribution, also caused by 
the spin-orbit interaction, is often referred to as the ex- 
trinsic spin-orbit coupling. In the first order with re- 
spect to the extrinsic spin-orbit coupling, there appear 
spin currents which are not equal to zero in the static 
limit. This leads to the extrinsic spin-HaU effect,^°'^^'^^ 
which is beyond the scope of the present paper. The ex- 
trinsic spin-orbit coupling also leads to an additional in- 
duced spin polarization, which is considered, in the limit 
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w ^ 3> rip, in Ref. 53. In the second order with re- 
spect to the extrinsic spin-orbit couphng there appear ad- 
ditional relaxation terms in the kinetic equation. These 
terms are responsible for the Elliot- Yafet spin relaxation 
of 2D electrons. ^^'^^ The corresponding spin relaxation 
rate is many orders of magnitude smaller that the mo- 
mentum relaxation rate and often can be neglected. 

The other important approximation is the spatial ho- 
mogeneity of the problem, which implies the neglect 
of the gradient terms in the kinetic equation (2). For 
this reason, the results of this paper cannot be di- 
rectly applied for description of the spatial distribution of 
spin density in 100 /um wide 2D layers recently studied 
experimentally.^'^^ Prom the theoretical point of view, 
the problem of finite-size samples is difficult, not only 
because of the presence of gradient terms, but also be- 
cause of the need to derive the boundary conditions for 
the matrix Wigner distribution function. In the general 
case, the required boundary condition should be a matrix 
equation."'''* Some simple forms of boundary conditions 
used in theoretical description of the spin-Hall effect (see, 
for example, Ref. 16) have been written without a deriva- 
tion. Therefore, the problem of spatially inhomogeneous 
distribution of the spin density still remains topical. 

The main qualitative result of this study is the pre- 
diction of vanishing spin polarization at arbitrary fre- 
quency of the applied electric field for the special situ- 
ations when fixed (momentum-independent) precession 
axes exist. The well-known examples of these situations 
are asymmetric [001]-grown quantum wells with equal 
Rashba and Dresselhaus velocities (vj^ = vjy) and sym- 
metric [110]-grown quantum wells {vr = 0). This re- 
markable property, which does not follow from the con- 
sideration of the static response, is explained by the ab- 
sence of DPK spin relaxation for these special situations. 
As a consequence, in the close vicinity of the fixed pre- 
cession axis, which means that v\ is close to for [001]- 
grown quantum wells or is close to zero for [110]-grown 
quantum wells, the frequency dispersion of spin polariza- 
tion is significant in the low-frequency region, when the 
polarization can be excited by applying an alternating 
electrical bias to the 2D sample. For [110]-grown quan- 
tum wells this dispersion is given by Eq. (63). Since 
the Rashba velocity is variable by modifying the shape 
of the quantum well via a bias applied to an external 
(top) gate, these phenomena can be investigated exper- 
imentally. As concerns possible technological use, these 
phenomena are of interest for the purposes of frequency 
filtering and amplification in future spintronics devices. 

It should be noted that the predicted behavior in the 
vicinity of the fixed precession axes can be influenced by 
the extrinsic spin-orbit coupling. From this point of view, 
the [001]- and [110]-grown quantum wells appear to be in 
quite different positions. In [001]-grown quantum wells 
with v\ = vjj, where the fixed precession axis lies in 
the plane of motion, a coupling of 2;-polarized spins with 
in-plane polarized spins occurs when the spin-dependent 
contribution to the scattering potential is taken into ac- 



count. This means that the applied electric field can ex- 
cite the spin density perpendicular to the precession axis 
as a combined effect of extrinsic spin-Hall current exci- 
tation and precession. Also, as mentioned above, there 
appears the Elliot- Yafet relaxation of the in-plane spin 
density. As a result, finite spin polarizations will ap- 
pear for the directions both perpendicular and parallel 
to the precession axis, and the non-analyticity with re- 
spect to the order of limiting transitions vj^ v'jj and 
a; — » will be removed. Still, the frequency dispersion 
of the spin polarizability should persist down to the fre- 
quencies comparable with the Elliot- Yafet relaxation rate 
uey- Using the estimate vey jv ~ 10^'"' given in Ref. 42 
and taking into account that v ^ 10*^ s^*, one finds 
Vey ~ 10'' s~*, which means that the frequency dis- 
persion can be principally observed in the region above 
10 MHz. On the other hand, in [110]-grown symmetric 
quantum wells, where the fixed precession axis is per- 
pendicular to the plane of motion, the applied electric 
field cannot excite the in-plane polarized spins even in 
the presence of the extrinsic spin-orbit coupling. Exci- 
tation of the current of i;-polarized spins (the extrinsic 
spin-Hall current) does not lead to spin polarization in z 
direction far from the boundaries of the sample. Finally, 
the Elliot- Yafet relaxation of i;-polarized spin density is 
absent in [110]-grown symmetric quantum wells. In con- 
clusion, the extrinsic spin-orbit coupling does not lead to 
a finite spin polarization in [110]-grown symmetric quan- 
tum wells and cannot modify Eq. (63). Modification 
of the low-frequency behavior described by this equation 
may occur owing to non-Markovian memory effects"''^ and 
effects of spatial inhomogeneity due to finite sample size. 

The other important result is Eq. (66), which es- 
tablishes a simple relation between the electric current 
and induced spin polarization. Though the validity of 
this equation is restricted by the approximation of short- 
range scattering potential, Eq. (66) is applicable to the 
whole class of the systems described by the p-linear spin- 
orbit Hamiltonians and does not require the linear re- 
sponse approximation. It can be useful in applications 
and in analysis of experiments where both spin polariza- 
tion and electric conductivity are measured. Ecjuation 
(66) also shows that the spin-orbit term in the electric 
conductivity is of purely dynamic origin, this term van- 
ishes at w ^ 0. It is not clear whether this property 
is specific for the short-range scattering potential or re- 
mains valid for arbitrary scattering potential. The cor- 
responding calculations are now in progress. In combi- 
nation with Eq. (43), the result (66) relates the elec- 
tric current with the spin current and leads to a unified 
description of spin and charge response to the applied 
electric field. 

The results for frequency dispersion of the spin po- 
larization and spin current obtained in this paper can 
be directly reformulated for description of transient spin 
response, ^^'^^'^^ which is investigated experimentally by 
the time-resolved spectroscopy.^'**''"'^ If the spin-split 
electron states are well-defined, the transient process 
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shows coherent oscillations. If the spin splitting is Acknowledgments 
suppressed by coUisional broadening, there should ap- 
pear long-time transients associated with the DPK spin 
relaxation. ^^'^^ Under conditions close to the appearance 
of the fixed precession axis, for example, when vr is close 
to zero in [110]-grown quantum wells, the duration of the 

transient process, according to Eq. (63), is expected to The author is grateful to Sergey Tarasenko for a helpful 
increase dramatically. discussion. 
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